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FIELD OF THE INVENTION 

[0001] Seismic traces of 3D surveys are interpolated using polynomials in 4D space 

from a less densely sampled real data set of source points S(x,y) and receiver points R(x,y). 
Input seismic traces should be pre-processed by nomial move out (NMO) and static 
corrections. Input traces contributing to each single interpolated trace are chosen from a 
selected number of traces closest to each interpolated trace in 4D space. The methodology 
in accordance with the invention is designed to minimize computational load. The amplitude 
of an interpolated trace is modeled with a two-temn Amplitude vs. Ofteet (AVO) fomnulation, 

BACKGROUND OF THE INVENTION 

[0002] In the oil exploration industry, seismic data is obtained to enable scientists and 

engineers to develop a picture of underlying rock formations. The reflection seismic method 
attempts to image the top few kilometres of the earth's crust by artificially creating a 
wavefield at the earth's surface and then recording this wavefield at multiple locations as it 
returns to the surface mainly via reflections within the rock of the earth's crust. These 
wavefields are then processed in order to obtain images of the subsurface that can be used 
to help locate hydrocarbons. 

[0003] As illustrated In Figure 1, a wavefield is created at the surface at a source 

location, S. by setting off a percussive shock wave or shot that imparts wave energy into the 
ground. The shot is typically dynamite or a mechanical vibration system, such as a Vibroseis 
truck, that creates a sinusoidal signal of changing frequency through shaking the earth. The 
energy source approximates a point source. 

[0004] A series of receivers or geophones at receiver points, R, set up as either a 

linear array (termed 2D seismic) or a two-dimensional array (3D seismic) around S. record 
the amplitude of wave energy that is received at each receiver point from reflections of the 
shot energy off underlying formations as a function of time, thus creating an array of 
time/amplitude data for each geophone location. In Figure 1 , schematically illustrating 3-D 

seismic, a two-dimensional array of receiver points ((1.1). (1.2). (1,3) ) is shown around 

the shot, S with representative wave paths that might be received at (1,1) based on 
reflections off the underlying fomiations. The recording of the wavefield is not continuous but 
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rather is sampled at discrete time intervals (typically 2ms intervals) and space intervals, 
which is determined by the separation distance between geophones (typically 50 meters). 

[0005] Reflections occur when the wavefield encounters a change in acoustic 

impedance, usually found at the boundaries between different rock types. Through various 
mathematical and statistical imaging techniques, an image of the sub-surface formations can 
be detennined at varying levels of resolution. As shown schematically In Figure 1 A, a number 
of traces with a certain density of spatial distribution can be used to develop a picture of the 
depth and shape of an underiying discontinuity. However, as can be seen from Figure 1A , 
no data exists for location A. 

[0006] From an imaging perspective, it is prefen^ed that the wavefield sampling is 

intense both in tenms of time (ie smaller sampling interval) and space (higher density of shots 
and geophones). In practice, however, the sampling is usually more than required in time, t, 
and less than required in space (x,y). In addition, azimuthal sampling is usually less than 
adequate (that is source to receiver azimuths tend to be similar rather than evenly distributed 
over 360 degrees) and the distribution of offsets (source to receiver distance) is unbalanced 
and inadequate. 

[0007] In general, therefore, the relative density of the shot and/or receiver points in a 

given area will enable the creation of images with relatively higher or lower resolution. 
However, as the cost of obtaining an ideally sampled dataset would be prohibitively 
expensive, the collection of seismic data must always be a balance between acceptable 
resolution and the cost of obtaining data with an acceptable resolution. Moreover, the 
collection of seismic data Is often affected by other limitations including the specific surface 
geography of an area where because of specific natural geographical features such as lakes 
or steep terrain, man-made features, such as roads or buildings, or environmental 
restrictions, such as wildlife sanctuaries, the placement of shots and geophones in particular 
areas is prevented, thus leading to incomplete images or Images with decreased resolution 
in particular areas. 

[0008] Accordingly, as the cost of signal processing is relatively inexpensive 

compared to the cost of field data acquisition, there continues to be a need to improve the 
resolution of the images obtained from seismic data using the industry acceptable densities 
of geophone arrays. In addition, there Is a need for systems that enable the interpolation of 
data into areas where a lack of shots and/or geophones have resulted in a lack of data. 



[0009] Other ProbI ms with Inadequately Sampled Data 

[0010] still further, a common manifestation of inadequate sampling is a low signal- 

to-noise ratio In the final image. Noise is a direct hindrance to the successful Interpretation of 
seismic data, and great efforts are made to reduce noise. One very successful way to 
eliminate random or quasi-random noise is to record redundant data. Using this method, all 
raypaths that reflect from a common point in the subsurface are added together or stacked to 
create a stadced trace of high fold. As the desired reflection energy adds constructively from 
each trace, and the random noise does not, the random noise is cancelled out. As a result, 
the higher the fold, the higher the signal-to-noise ratio. 

[0011] As noted above, in order to con-ectly image a geological structure in the 

subsurface, the surface spatial sampling must be adequate. Sampling theory describes or 
explains that the steeper an underlying structure is, the tighter the spatial sampling on the 
surface must be in order to correctly image It If the spatial sampling is inadequate (le 
geophones are too far apart), then spatial aliasing will result giving Incorrect images. 

[0012] Similarly, imaging complex structures may also require adequate sampling 

from a variety of azimuths to ensure that a variety of raypaths as defined by source-to- 
receiver raypaths provide sufficient data to avoid the reflection problems realized by complex 
structures that would otherwiseresult in a distorted or incomplete image. 

[0013] Further still, the reflection amplitude and phase are not constant at all angles 

of reflection where both reflection amplitude and phase will vary with the angle of inddenoe 
of the raypath. These variations provide important Information about the rock layers at an 
interface. In addition, various forms of unwanted coherent noise (as opposed to random 
noise) travel at velodties that differ from the primary (and desired) reflection energy. This 
difference in recorded amval times becomes larger with increasing distance traveled (as 
defined by source-to-receiver offset) and thus can be used to effectively remove the 
unwanted coherent signal. However, in order to do this, a good balance of source-to-receiver 
offsets must be present at all points in the recorded dataset. 

[0014] Further still, various processing steps rely on good statistics throughout the 

dataset. A less than adequate dataset from the field can limit the quality of the image 
processing. Steps that are hampered by inadequate sampling include derivation of refraction 
statics, deconvolution. operator design, reflection statics detemiination. velocity analysis, 
various noise attenuation techniques, pre- and post-stack migration and alt reservoir analysis 



techniques, such as AVO/LMR analysis that use the variation of reflection amplitude with 
offset as their input. Shortcomings introduced by inadequate statistics in any of these steps 
can seriously impact the final image integrity. 

[0015] As a result, there is a need for a system that continues to address the 

problems identified above with respect to the collection and interpretation of seismic data. In 
particular, there has been a need for a system that conveniently and effectively enables the 
interpolation of seismic data to provide an accurate or improved resolution of underiying 
formations around which there may be a shortage of data. 

SUMMARY OF THE INVENTION 

[0016] The invention seeks to assist in predicting the fomriation structure at a location 

where a shot or geophone may not be present by providing a method for interpolating data 
obtained from adjacent shots and geophones. 

[0017] Generally, the technique involves fitting a multi-order polynomial equation to 

the data from multiple locations around the desired location. 

[0018] More specifically, there is provided a method of generating output seismic 

trace data at a specified location from input seismic trace data from known locations 
comprising the steps of selecting a spedfied location for the output data, selecting input data 
at known locations for fitting to a multi-order polynomial equation, obtaining polynomial fit 
parameters of the multi-order polynomial equation and utilizing the polynomial fit parameters 
to obtain output data at the selected location. 

[0019] In further embodiments, the step of selecting the location for the output data 

includes selection of both a desired source point and a desired receiver point. 

[0020] In another embodiment, the step of selecting input data at known locations for 

fitting to a multi-order polynomial equation includes selecting input data within a specified 
range of the specified location of the output data. 

[0021] In a still further embodiment, the input data within the specified range is further 

restricted by selection of input data closest to the specified location of the output data. 

[0022] In one embodiment, the step of obtaining polynomial fit parameters of the 

multi-order polynomial equation is solved numerically by a singular value decomposition of a 
matrix, C, defined by C = ULV^ and in another embodiment a transformed coordinate 
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system Is utilized to speed up the computation time required to solve for the polynomial fit 
parameters of an interpolated trace. 

[0023] In yet another embodiment, each of the steps is repeated to fomn a shot 

gather corresponding to predicted receiver trace data at known or unknown locations and 
corresponding to a known source location or predicted source location. 

[0024] In another embodiment, the specified location is a known location and the 

method is utilized as a means of reducing random noise associated with the trace data at 
that known location. 

[0025] In another embodiment, the invention provides a method of creating a shot 

gather from known seismic trace data obtained from a plurality of known locations within an 
x,y grid, the method comprising the steps of: 

selecting a source location within the x,y grid to which the shot gather will 
correspond; 

selecting a plurality of receiver positions that will fomi receiver positions for 
the shot gather wherein at least one of the source location or receiver 
positions is an unknown position; 

selecting a plurality of receiver traces having known positions as input traces 
and utilizing the known receiver traces to predict receiver trace data at a 
location adjacent the known receiver positions utilizing a least-squares fitting 
model; 

repeating the third step for all selected receiver positions so as to fomn a shot 
gather corresponding to the source location. 

[0026] In a further embodiment, the source location is a location on the x,y grid that 

does not correspond to an actual source location. 

BRiEF DESCRIPTION OF THE DRAWINGS 

[0027] The invention is described with reference to the following drawings wherein: 

[0028] Figure 1 is a schematic isometric diagram of an array of receiver points, R, 

around a shot point, 8, showing possible raypaths to receiver point (1,1). 



[0029] Figure 1A is a schematic diagram of part of a regular receiver grid at fpur 

receiver points; 

[0030] Figure 2 is a flow chart of the main steps in accordance with the invention; 

[0031] Figure 3 is a schematic diagram of input shot positions on a 96 by 513m grid; 

[0032] Figure 4 is a schematic diagram of a line receiver patch having 7 receiver 

lines and 128 receivers; 

[0033] Figure 5 is a schematic diagram of receiver patch size of 128 by 13 with 

added stations; 

[0034] Figure 6 is a schematic diagram of shot grid size of 56 by 21 with added 

stations. 

[0035] Figure 7 a surface map detailing the location of real shot and receiver lines; 

[0036] Figure 8 shows a representative shot gather for real shot 388 for every third 

live receiver station. 

[0037] Figure 9 shows a predicted shot gather for real shot 388 wherein 50% of the 

real shot gather traces were assumed not to exist; 

[0038] Figure 10 shows the difference between the predicted data of Figure 9 and 

the raw data of Figure 8 for shot 388 and represents the lack of prediction of random noise; 

[0039] Figure 1 1 shows a predicted shot gather for shot 388 In which the system was 

used to predict the shot gather utilizing all real shot gather traces; and, 

[0040] Figure 12 shows the difference between the predicted data of Figure 11 and 

the raw data of Figure 7. 



DETAILED DESCRIPTION OF THE INVENTION 

[0041] With reference to the Figures, a system for the interpolation of seismic data 

utilizing known seismic data is described. The system provides a method of processing 
known seismic data for predicting seismic data at locations for which no data is known 
utilizing a multi-order polynomial fit methodology. 

[0042] In accordance with the Invention, seismic data (amplitude, time) obtained from 

an x,y grid of known seismic data is utilized to obtain predicted data arrays at locations within 
the grid where no data has been obtained. As shown in Figure 1 A, the invention enables the 
prediction of seismic data at a receiver location A for which no data has been collected. The 
methodology of the invention generally includes the steps of preparing the input data 10. 
selecting the location of the output or predicted data 12. selecting the location of the nearest 
location input data 14, fitting the input data to a multi-order polynomial equation 16 and 
thereafter utilizing the polynomial parameters of the fit equations 18 to obtain predicted data 
arrays at the selected locations. The basic processing steps are repeated in order to prepare 
shot gathers that present trace data corresponding to known or unknown shot locations and 
known or unknown receiver locations. The methodology is also effective in providing shot 
gather data that has been filtered to eliminate random noise associated with actual trace 
data. 

[0043] Preparing Input Data 10 

[0044] Initially, it is preferred that unstacked input data Is pre-processed In 

accordance with known prestack signal processing steps including but not limited to 
deconvolution, amplitude scaling* statics, velocities, noise attenuation, spherical divergence 
corrections, static conrections and normal moveout (NMO) corrections and others as may be 
known to those skilled in the art. Pre-processing may take place at any stage prior to fitting 
16. 

[0045] Select Output Data Locations 12 

[0046] The location of a new predicted data point (one single time sample) at a new 

output location Xs, ys» Xr, yr is selected. The selection of the location of the predicted data 
point may be based on a desire to increase data resolution between lines of actual data or by 
a desire to predict data at a location from which neither a shot or receiver could be 
established by virtue of geographical limitations. 



[0047] More specifically, it is preferred that selecting the output data location involves 

defining a) the desired source point and b) desired receiver patch locations where a 
prediction is desired and defining the x,y regions where you have both c) actual source point 
and d) actual receiver point data. The desired source point and receiver patch locations are 
preferably in the form of a rectangular grid. 

[0048] Select Input Data Locations 14 

[0049] All input data amplitudes within a specified distance of the new output's x.y 

location and at the same time are defined as input samples. The input data selection will 
preferably be within a formula based range such as an ellipse or circle, although other 
ranges may be defined. 

[0050] In greater detail, and in a preferred embodiment, input data is chosen with 

"nearby" common mid points (CMPs). In this embodiment, for every output data trace with 
source location at (Xs.ys) and receiver location at (Xr,yr), input traces are chosen for 
interpolation. An input trace with source location at (Xs.ys) and receiver location at (Xr,yr) may 
satisfy a specific formula based condition such as that which defines all input data within an 
ellipse (or other shape such as circle), 

[0051] (x^-^^)' ^ (y« Y'^^' ^1 (1) 

[0052] where Xm=( Xs+ Xr)/2. ym=( ys+ yr)/2. Xm=( Xs+ Xr)/2. y^={ ys+ yr)/2. and r^ and Ty 

are user specified sizes, the larger the more spatial mixing effects will occur. 

[0053] The data may be further restricted to choose data in "nearby" 4D space. For 

example, from the above, there may be many input traces satisfying the proximity conditions 
defined by an ellipse. However, in one embodiment, only k input traces will be kept, with k 
being defined as the kept trace number or user specified input fold. The higher the input fold, 
the cleaner the interpolated traces will be. The kept traces (k) are the ones the closest to the 
desired trace to be interpolated in 4D space, defined by distance in 4D. 

[0054] -^,)' -*-(y, -x,r +(y, -y^rr' 

= rf4(x,,y^,x^,y^,jr, ,x^,y^) 
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[0055] Only the k input traces with the smallest d4 will be kept with the input traces 

being ranked for selection. The total of k input traces have their source & receiver location 
denoted by 

[0056] x^U\ysUUrU).yrU\ J = iA..Jt (3) 

[0057] Practically, k=15-100 input traces with 15 being the minimum to be able to 

solve a 2"^ order polynomial In 4D space. The selection of k is a trade-off between 
computation time and the benefits of obtaining a smoother Interpolated trace that Is less 
corrupted by noise and mixing effects. 

[0058] Fit Multi-Order Poiynomiai 16 and Predict Data for Seiected Location 18 

[0059] A multi-order polynomial is fit to the amplitudes of the input samples. 

Practically, a second-order polynomial of 15 coefficients is usually adequate. However, a 
polynomial of higher order is generally prefen^ed for more complex structures, steeper dips, 
when the 4D distance between Input traces and desired output location becomes large and 
to accommodate complex AVO effects. 

[0060] At a time sample t. t=dt, 2dt, ... Ndt. the jth input trace has amplitude. 

A{t.x,ijXy,(AxAjXyrO}}^Ait,Xj) (4) 

[0061] It Is fitted to a polynomial in the 4D coordinates: 

[0062] />(%)= Z^y^^i^/^rJ^'i ^5) 

[0063] where, m & M are user specified single and total power numbers, the higher m 

& M, the more accurate the input amplitudes can be fitted, and the higher the resolution of 
the fitting. The polynomial coefficients, a^f^ are solved by minimizing: 

[0064] = - P{xj)r . (6) 

[0065] This is equivalent to solving the least squares solution to the set of equations: 

[0066] P(Xj ) = A(t, Xj ), J = 1,2,.../: (7) 



[0067] The set of equations above are linear in unknown polynomial coefficients a^^ . 

With a^ia solved, interpolation is perfonmed by substituting the coordinates of the trace to be 
interpolated into the polynomial: 
[0068] P(x) = P(x„y„x,,y,). (8) 

[0069] The amplitude of the new output location is predicted by the polynomial at the 

new output location x,y coordinates. 

[0070] The time sample at t of the pre-stack trace at (Xs,y8.Xr,yr) has now been 

interpolated. 

[0071] The steps are repeated for all time samples in order that the pre-stack trace at 

(X6>y8>xr.yr) is interpolated and complete arrays of interpolated shots are created. This 
completes the interpolation of the full data set. 

[0072] The system is preferably implemented in an output-oriented way enabling 

user-selection of data-selection and data-handling parameters as explained in greater detail 
below. 

[0073] System Optimization 

[0074] In various embodiments, the system may be optimized in accordance with the 

following methodology. 

[0075] Preparation 

[0076] Write equation (7) in matrix form: 

I 



[0077] 



I 
I 



[0078] or, 

[0079] Ca = A (10) 
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[0080] In (9). each equation represents one input trace having an arbitrary order. In 

the vector a, the order of polynomial coefficients is also arbitrary, but for convenience, its 1^ 
element is aoooo- 

[0081] The matrix C has its singular value decomposition: 

[0082] C = ULV^. (11) 

[0083] where superscript T denotes matrix transpose, and L is a diagonal matrix of 

positive singular values, (...U...). The singular value decomposition of C can be easily 
solved by many readily available numerical software packages. 

[0084] The generalized Inverse of matrix C, denoted by C"* . is given by 

[0085] C-«=VL-'U^. (12) 

[0086] where, L*^ is a diagonal matrix of (...l/IU+s^. •). and is a small positive 

damping, the effect of which is to stabilize the matrix inversion and at the same time 
interpolated results will be smoother and cleaner. This generalized inverse provides a 
solution to polynomial coefficient vector, 

[0087] a = C~*A (13) 

[0088] This is a least squares solution when matrix C has more rows than columns. 

[0089] To speed up the solution of equation (9) or (7). since matrix C is not time 

dependent, is also time independent. C*^ needs to be calculated only once for all the time 
samples of this pre-stack trace to be interpolated. 

[0090] This reduces computation to 1/(number of time samples per trace), typically 

1000+ time samples. 

[0091] To speed up the calculation of equation (8), the origin of 4D space is moved to 

the position of the trace to be Interpolated, i.e., replace {x^,y^,x,,y,) by (X8-X8,y8-ys,xrXr,yryr). In 
this transformed coordinate system, the position of the trace to be interpolated becomes 
(0,0,0,0). Interpolating the trace by substituting its position into the polynomial, only the zero 
power term survives. In other words, the interpolated sample is the constant term of the 
polynomial, i.e. aoooo- Therefore, instead of solving the full vector of polynomial coefficients, 
and calculating the value of the polynomial at the location of the trace to be interpolated, we 



only need to solve aoooo, in the transfonned coordinate system. From (13) it is seen that, only 
the row of corresponding to aoooo needs to be calculated, and 

[0092] aoooo = C-9(1 ^ row) * A , 

[0093] where, denotes dot product of two vectors. 

[0094] This reduces computation of polynomial coefficients to 1/(number of 

polynomial coefficients), and completely removes the need of calculating the value of 
polynomial at the output position. 

[0095] To speed up the whole process of interpolating one pre-stack trace at Xs, ys. 

Xr. yr for all t-samples. combining the above two, the interpolated trace is a weighted sum of 
the input traces with the weights being the elements of the 1st row (comesponding to aoooo) in 
the time independent C^*: 

k 

[0096] ^(^x,,y,,x^,y^) ^^C;fA(t, x.uhysUhX^uhyrUi) , (14) 

[0097] in which C^f is time Independent and is a function of input trace positions (x, 

/) and the position of the interpolated trace (x, y). 
[0098] User Options 

[0099] As indicated above in various embodiments, it is preferred that the user be 

able to specify various options including: 

[00100] 1. Output Shot Gathers 

[00101] A rectangular shot grid in (x,y) including shot start/end x and y coordinates 
and how many shots between start/end x and y coordinates as well as a rectangular receiver 
patch grid relative to each shot of the shot grid. 

[00102] 2. Output Common Mid-Point (CMP) Gathers 

[00103] A CMP grid In (x, y) including CMP start/end x and y coordinates and how 
many CMPs between the CMP start/end x and y coordinates. In addition, it is preferred that n 
offset/azimuth patch for each CMP gather is enabled including a shot/receiver offset 
(minimum and maximum and how many in equal interval) and shot/receiver azimuth 
(minimum and maximum and how many in equal interval). 
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[00104] In both cases, in order to avoid extrapolation, not all points in the grids can 
either have a live shot/receiver or within some CMP gathers it may be impossible to create 
some or all of the desired offsets/azimuths. 

[00105] Examples 

[00106] A number of tests were performed on both real and synthetic data. 
[001 07] Synthetic data 

[00108] Testing using synthetic data utilized a delta function approach. A zero phase 
wavelet is reflected from a tilted boundary with random source and receiver locations about a 
common mid-point (CMP). The gathered traces contain the same wavelet placed at 
calculated arrival times. With such an input data set, interpolation to the same locations has 
resulted in near perfect results. 

[00109] Field Data 

[00110] A few field-data sets have been utilized that provide satisfactory interpolation 
results as detailed below. 

[00111] Sources fonm a grid with x and y intervals of dx=513 & dy=96 metres, with 6 
shot lines oriented in the y direction, and 56 shots per line. Note that geographical features 
and dead shots either prevented a shot from being made or resulted in a no shot. The total 
number of live shots Is 317 as shown in Figure 3. 

[00112] The shot data was received with receivers fomning a grid with x and y intervals 
of dx=16 & dy=480 metres, with 12 receiver lines oriented in the x direction, and 160 
receivers per line (not shown). 

[00113] A live receiver patch for each shot has 7 receiver lines and 128 receivers per 
line as shown in Figure 4. Thus, there are 896 traces per input shot. 

[001 14] Both source and receivers occupy the same rectangular zone of 2550 by 5280 
(x by y) metres. 

[00115] In accordance with the invention, interpolation tests were perfomied to: 
[001 16] a) Expand the receiver patch/grid 

[00117] As shown schematically in Figure 5, the number of receivers per line (128) 
was k pt the same, but an additional receiver line was inserted, via interpolation, between 
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the existing 7 receiver lines of the patch thus increasing the number of receiver lines in the 
patch to 13. This is equivalent to reducing the receiver line spacing from 480 to 240 metres. 
All interpolated receivers were within the rectangular zone of existing receivers. 

[00118] b) Expand the source grid 

[00119] The source grid was also changed by inserting through interpolation an 
additional 3 shots line between the existing 6 shot lines as shown in Figure 6, thus increasing 
the number of shot lines to 21. This is equivalent to reducing the shot line spacing from 513 
to 128 metres. 

[00120] c) Create Interpolated shot gather 

[00121] Utilizing the interpolation system in accordance with the system, shot gathers 
were created for each source grid shot (both actual and interpolated) having individual 
receiver traces for all receiver locations. 

[00122] d) Results 

[00123] Visual inspection of the interpolated shot gathers shows good results. Most 
interpolated data are stacked and/or pre-stack migrated. Visual inspection of the stacked 
sections and/or time slices shows preserved structural details, with increased S/N ratio, 
reduced acquisition footprints and a somewhat spatially smoothed look, which can be largely 
controlled by the fold parameter (k) and the damping (e^ ). 

[00124] Examples cyf Output Shot Gathers 

[00125] With reference to Figures 7-12, examples of shot interpolated and non- 
interpolated shot gathers are described In addition to other processing techniques to show 
the effectiveness of the method. 

[00126] Figure 7 shows a surface map detailing the location of real shot and receiver 
lines. In this example, the shot and receiver lines were placed along existing geographical 
features such as roads and rivers and, hence, the shot/receiver line density is substantially 
less than many seismic surveys. Thus, this example represents a relatively crude survey. 

[00127] Figure 8 shows a representative shot gather for real shot 388 for each real 
receiver station and, accordingly shows a typical shot gather as is known in the art. 

[00128] Figure 9 shows a predicted shot gather for real shot 388 wherein 50% of the 

real shot gather traces are assumed not to exist. Specifically, in generating the shot gath r, 
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each even numbered shot gather was ignored and the system in accordance with the 
invention utilized to predict the new shot gather traces at those even numbered location. As 
can be seen, in comparison to the raw data of Figure 8, the system is effective in predicting 
the true formation and does not model the random noise which is apparent in Figure 8. 

[00129] Figure 10 shows the difference between the predicted data of Figure 9 and the 
raw data of Figure 8 for shot 388. This shot gather represents the lack of prediction of 
random noise. 

[00130] Figure 1 1 shows a predicted shot gather for shot 388 in which the system was 
used to predict the shot gather utilizing all shots. As can be seen from Figure 12. which 
shows the difference between the predicted data of Figure 1 1 and the raw data of Figure 7, 
the use of the system to predict a shot gather at non-interpolated positions is effective in 
eliminating the random noise associated with actual data. 

[00131] Other 

[00132] When survey data is collected sparsely. Irregularly and with "holes", pre-stack 
Interpolation may be used to resample the data onto a regular and/or denser shot and/or 
receiver grid. The resulting benefits include: 

[00133] a) Stacking: Stacking fold is higher, leading to cleaner stacked sections. In 
addition, offset distributions are more unifomi, leading to less acquisition foot prints and more 
reliable offset range limited partial stack. 

[00134] b) Pre-stack migration: There are fewer dead traces, leading to less migration 
"smiles". In addition, the more finely sampled data makes numerical integration (the essence 
of migration) more accurate. 
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